B — > p form factors including higher twist contributions and reliability of pQCD 
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We discuss B — > p form factors within the framework of perturbative QCD, including the higher 
twist contributions and study the validity of such an approach in calculating quantities such as 
form factors, which in principle and quite generally are thought to be completely non-perturbative 
objects and are expected to receive large contributions from the non-perturbative regime in the 
calculations. It is shown that including the Sudakov and threshold resummation effects, the general 
expectations of the pQCD approach are clearly met and the form factors do indeed receive most of 
the contribution from the perturbative region. We do not make an attempt to precisely evaluate the 
form factors but rather try to study the gross features and behaviour of the same. We also find that 
use of single wave function for the B-meson may actually underestimate various quantities. The 
results clearly indicate the validity and reliability of pQCD calculations, at least in this particular 
case. 

PACS numbers: 13.20.He, 12.38.Bx, 12.38.Cy 



I. INTRODUCTION 

B-decays offer a very fertile soil for understanding, analysing and testing the basic structure of particle interactions, 
widely ranging from electroweak aspects, including CP violation, to the mysterious world of QCD (for a quick 
review of various issues see Q). The study of B-decays and related observables has greatly enhanced and affected 
our understanding of the underlying principles that guide and govern these phenomena. The Standard Model (SM) 
of particle physics, including the QCD corrections, seems to be a highly successful candidate in explaining almost 
all the experimental data, including a variety of issues concerning B-decays themselves. However, as we enter the 
precision era in B-physics, there is a compelling need to uncover more and more that goes as input while evaluating 
and analysing these decays, and as can be expected, the need to unearth the relative importance of QCD corrections 
is overwhelming. In particular, the experience with b — > 57 shows that QCD can, in fact, alter the results by a large 
amount and therefore it becomes imperative to include higher order QCD corrections to get more and more sensible 
and accurate results. 

Scmileptonic decays of the B-mesons are of particular interest owing to their cleanliness and relatively simpler 
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calculational aspects. It also means that these decays provide us with the opportunity to test QCD corrections 
more reliably and much more easily, compared to pure hadronic decay modes. This has been widely recognized and 
a lot of work has been done in this direction. Owing to the large mass of B-mesons, it is expected that the heavy 
quark effective theory (HQET) description is a good one and the leading term is given simply by the quark level 
process (see Q). The sub- leading terms are suppressed by the B-meson mass and in most of the applications, are 
therefore neglected. Another simplifying assumption that is employed is the idea of factorization |3|. Quite simply, 
it means that in an energetic process, like the decay of a B-meson, since the energy released is large, the outgoing 
hadrons move out of the interaction region quite quickly, thereby implying that there are no soft interactions 
between the various hadronic subsystems. The full meson level amplitude for any process is thus written as a 
convolution integral of the hard scattering kernel (which can be reliably computed using standard perturbative 
formalism) and the non-perturbative meson wave functions that are universal and are obtained from lattice studies, 
fits to the experimental data or sum rules of one kind or the other. This assumption forms the backbone of most 
of the calculations involving B-mesons. The issue of the extent of reliability of such an approximation, which seems 
to work quite well for most of the cases, is still an open question. It is however found and argued that not in all 
cases this is a good approximation and one must look for the sub-leading terms to this approximation scheme. For 
example it suffers from the problem of scale dependence jj] and in case of B — > J/^> it seems difficult to match the 
branching ratios using the approximation Q Furthermore, it is now accepted that though the leading term in the 
HQET gives results that are quite close to the observed numbers, the reality that the mass of a B-meson is not 
infinite but close to 5 GeV, implies that the sub-leading terms can not be naively neglected anymore. One can try 
and estimate these sub-leading terms in HQET itself or employ some other methods. 

It is worthwhile to try and explore the idea of estimating the above mentioned sub-leading terms suppressed by 
the B-meson mass. To do such a calculation, one can follow either of the following two approaches. Both the 
approaches treat the hard scattering kernels perturbatively and the mesonic wave functions are the non-perturbative 
ingredients. The main difference lies in the treatment of quantities like from factors. 



(a) QCD factorization: The form factors are also thought to be purely non-perturbative objects, not calculable 
perturbatively 

(b) pQCD: It is believed that the form factors can be reliably and satisfactorily computed in perturbation theory. 
This approach was developed by Brodsky and Lepage and others 7\. The factorization theorem for relevant 
for exclusive B-meson has been proved [g. 

There is no clear consensus on either of these two methods and both seem to have some advantageous features while 
both suffer from lack of very sound theoretical footing on certain issues. Both the approaches have invited intense 
activity in the recent past. In some cases, the results of the two approaches tend to be consistent while in some other 
cases, they seem to be distance apart. 

In this note, we investigate B — ► p form factors in the second of the approaches listed above, namely the pQCD 
method. We have chosen semileptonic process for the merit of cleanliness and relatively simplified calculations and 
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the fact that with just two hadrons in the process, we do not have to bother about the subtleties and complications 
involved in a purely hadronic process and the reliability can be checked more clearly. The same process has been 
considered by 9] and io| . We extend their study and elaborate on the differences in the next sections. The aim 
of this study is not to precisely pin down the value of the form factors from the calculations but to explore and 
study the behaviour of the same with respect to the various parameters that enter the calculation and check the 
reliability of such a calculation. In particular, we keep in mind the objections and criticisms that are generally raised 
regarding the validity of such a scheme and try to see for ourselves, whether, following a more consistent treatment 
compared to whatever exists in literature till now, we get a clue to some of the unresolved and mysterious issues 
that we are forced to live with in such calculations. We would like to emphasize again that the goal is not a very 
accurate numerical study of the form factors but to study the general behaviour for some suitable choices of various 
parameters like the shape variables appearing in the meson wave functions etc. Furthermore, we would like to remind 
ourselves that in such a situation, we may finally end up over- or under-estimating some of the quantities but we 
hope to clarify certain issues in the end. In this spirit, this study aims at extending the ongoing debate between the 
two approaches in order to get a clearer picture of what exactly is happening and can we understand and explain 
the same. The article is organized as follows: in the next section we briefly review the chief ingredients of the pQCD 
approach. We summarize some of the main objections/criticisms that this method faces. Next we discuss the form 
factor calculation including the higher twist contributions to the meson wave functions. We then qualitatively study 
the behaviour of the various form factors without attempting to determine the values very precisely. The last section 
discusses and summarizes our results and conclusions. 



Exclusive processes enjoy the simplicity and elegance they derive from factorization theorems invoked in some form 
or the other (for a general review of factorization theorem and its related issues see 11]). It has been demonstrated 
that for an exclusive process involving a large scale Q, the amplitude can be neatly written as 



where C(t) denotes the Wilson coefficients relevant to the problem, H{t) is the hard scattering kernel that is per- 
turbatively evaluated, describes the distribution of the partons in the i-th hadron (here x is the momentum 

fraction carried by the parton) and we have collectively put all the resummed quantities in the factor S for simplic- 
ity. The factor S therefore contains the Sudakov logarithms and the relevant evolution factor (these are discussed 
below). The factorization theorem implies that in a hard exclusive process, the non-perturbative dynamics can be 
separated from the perturbative pieces and the final result is simply the convolution of these. We now concentrate 
on a specific process involving B-meson decay, namely B — > ptv, as our prototype hard exclusive process. The large 
mass of the B-meson acts as the natural large scale in the problem. Here and in the following we do not differentiate 
between a B or a £?-meson. We therefore have the following scales in the problem: mass of the b-quark (we do not 
differentiate between the b-quark mass and the B-meson mass) rr% ~ 5 GeV, the W-boson mass Mw ~ 100 GeV, 
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(2.1) 
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the renormalization scale p (which is of the order of nib) and a scale of the order of Aqcd = A ~ 250 GeV. The 
Wilson coefficients appearing in the effective Hamiltonian contain the information of the high scales (> nib) and 
have the resummed logarithms of the form ln(Af^//i 2 ). The scale A characterizes the non-perturbative scale and 
crucially enters the meson wave functions. The precise shape and behaviour of the meson wave function are crucial 
inputs for any calculation to make physical sense. For example, if the wave function does not vanish at the end 
points (x — > 0, 1), then it can be shown that the amplitude is infra-red divergent [l^ . Also, it has been pointed out 
that if the small transverse momentum of the parton, generically denoted as k± ~ 0(A) is ignored, the dominant 
contribution arises from the end point regions. The way out is to retain the k± components, which in turn regulate 
the infra-red divergence that has just been mentioned [ljj- If the variable conjugate to k± is denoted by b, then 
apart from the large logarithms that are reorganised using standard renormalization group techniques into Wilson 
coefficients, we have another potential source of large logarithms - \n(pb). The inclusion of the transverse momentum 
for the partons in turn results in double logarithms of the form ln 2 (P6), where P is the momentum of the meson 
(generally the larger of the light-cone momentum components). These large logarithms are resummed and lead to 
the so called Sudakov form factor exp[—S(p,b)]. This exponentially damping factor suppresses the long-distance 
contributions from the large b regions and the factor S(P, b) vanishes as b — > 1/A. Since the infra-red divergences in 
a theory are manifestations of the non-perturbative dynamics, such effects are absorbed in the wave functions. These 
wave functions are universal in character and do not depend on the specific process. After, reorganising all the large 
logarithms and absorbing the infra-red divergences in the meson wave functions, we are left with finite quantities 
that are believed to be perturbatively calculable. The main philosophy behind the pQCD method is that all the 
factors except the hadronic wave functions can be perturbatively calculated such that the dominant contribution to 
various quantities, like form factors, comes from the hard gluon exchanges. It is important that this is true because 
otherwise the dominant contribution would arise from the non-perturbative regime - a region where the perturbative 
calculations do not make sense and there is absolutely no way of determining them in this way. 
The B — > p form factors are parameterized as follows [b^ : 

(p(P p ,X)\(V-AUB(P B )) = -iimB + mJAitfK™ + lM / ] ( e * ■ P B )(P B + P p ), (2.2) 

(TUB + nip) 

iA 3 (q 2 ) 2V(q 2 ) p „ pa *0 
+ 7 ; d e ' Pb)(Pb - P P )n + t ; ^"otpPB^p e ± 

where q 2 = (Pb — P P ) 2 is the momentum transferred. For the process under consideration, we have the following 
two diagrams contributing at the one gluon exchange level in pQCD: 
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FIG. 1: One gluon exchange diagrams contributing to B — > p form factors. The cross denotes the weak interaction vertex, 
where <pB,p represents the corresponding meson wave function, £ is the momentum fraction carried by the lighter 
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(anti)quark and l±_ is the transverse momentum while x and k± denote the analogous quantities for the p-meson. 
As mentioned before, we do not distinguish between b-quark mass and B-meson mass and set them equal in the 
calculations. The light quarks - u,d, s are taken to be massless. Our convention for the light-cone variables is as 
follows. For any four momentum = (a + , a~ , ajj we define a ± — a ± a 3 and a^ = (a 1 , a 2 ). The scalar product is 
defined as: a ■ b = a b ~ + a h - a ± b ± . With these conventions, we have for various momenta: 

P B =m B (l,l,0 1 ) P p = ^(r 2 ,7? 2 ,0 ± ) (2.3) 

V 

where r — m p /mB and rj = 1 — q 2 /m 2 B . It is convenient to define two light-like momenta n + = (V2, 0, 0j_) and 
n_ = (0, V2, Ox) such that n+ • n_ = 1. The transverse momenta k± and l± as introduced above should be thought 
of as four-vectors with only the transverse components non vanishing. Schematically, therefore, the various form 
factors (generically called F), after using the factorization theorem, can be written as 

F= [ dx [ dt<j> B {Z iq 2 )H{x,^q 2 )<p p {x,q 2 ) (2.4) 
Jo Jo 

where as before <^>'s denote the meson wave functions and H{x^,q 2 ) represents the hard kerne evaluated perturba- 
tively. The basic theme behind such a factorization is that the long-distance or the soft interactions occur before 
and after the hard decay process and therefore the two effects simply decouple. The hard contribution takes place 
at short-distance scales and therefore we only need to specify the distribution of quarks (partons) inside the mesons 
apart from explicitly computing the perturbative piece. The information regarding the parton distribution, along 
with all the infra-red dynamics, is contained in the wave functions. In the asymptotic region, the factorization 
theorem has been shown to be valid. However, it has been pointed out |12| that in the case of pion's electromagnetic 
form factor at lower energy scales (few GeV's), large contributions come from the end point regions (x,£ — > 0). The 
perturbative analysis is not valid in such a region and this casts serious doubts about such a method. This problem, 
even when the total convolution integral does not contain any divergences, is called the endpoint problem. The way 
out is to introduce Sudakov factors that regulate this undesirable endpoint behaviour and thus make the pQCD 
calculations sensible at these lower energy scales The Sudakov factors typically suppress the long-distance 

effects at large transverse distances (denoted by variable b earlier) ie. small transverse momenta. Therefore, the 
endpoint behaviour is a crucial aspect to get any meaningful results. 

A further issue of concern is that as the momentum fraction carried by the spectator tends to zero, the form factors 
become divergent due to divergent behaviour of the meson wave functions related to the light-cone distribution 
amplitudes (LCD As). There is a need to resum the contributions of the form ln 2 x. This is achieved by threshold 
resummation into a jet function St(x) such that St(x) — > as x — > 0, 1 16] . Therefore, the threshold resummation 
modifies the endpoint behaviour of the distribution amplitudes (DAs) . 

After including all these factors in the DAs, one expects that the pQCD calculations are reliable and can be carried 
out without much trouble. However, the method still faces some more objections/criticisms which we now briefly 
discuss. We summarize the main issues as discussed in \v\ 



Are Sudakov suppression factors strong enough to regulate the large transverse separation contributions? In 
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particular, it is quite possible that for intermediate values of the variables, away from the endpoints, there can 
be a slight enhancement. How effective and efficient is this effect? 

• How small are the contributions from the non-perturbative regime? 

• The biggest uncertainties presumably arise from the precise lack of knowledge about the meson wave functions. 
In particular, the use of a single wave function for the B-meson has been questioned. 

All these and related issues cast a serious cloud of doubt on the applicability and reliability of pQCD calculations. 
We try to address and clarify some of these issues below and also elaborate on various issues. 



The DAs refer to the distribution of partons inside a hadron (here mesons). A DA thus, carries the information 
regarding the momentum fraction carried by partons inside the meson in a specific Fock state. The DA is related to 
the the Bethe-Salpeter wave function by the following relation 



where /i is the ultra-violet cut-off (for a general discussion sec 18] . In the discussions below, we do not remain very 
careful and do not distinguish between the wave functions and DAs and freely interchange one for the other in the 
discussions. 

We closely follow the approach and results for the p- meson DAs as outlined in . Accordingly, the the LCDAs are 
defined as meson-to- vacuum transition matrix elements of the non-local gauge invariant operators in the light-cone 
picture. For the light vector mesons, the DAs are split into chiral even and chiral odd contributions. In terms of 
n±, the vector meson momentum is written as P p = En^ + m 2 p n + ^J (AE), such that the transverse plane is defined 
with respect to the vectors n±. Fourier transforming the DAs, we obtain the momentum space representation of the 
p-meson light-cone projection \2(\ (keeping terms to twist-3 of the two particle quark-antiquark distribution only 
and not writing the path ordered integral and also remembering a relative —i factor between the definitions here and 



III. MESON DISTRIBUTION AMPLITUDES 




(3.1) 



those in [19|) 



dc(*)t*«(0)|0) = M p s a = + M£ 
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M L± = if "'"" 



4 * v 4 



/^(x) - 25 / <M0||( V )-.g[%)] A-el"^ (3.4) 



\ ) J / 5a 

By <9 fei we mean and we have 

= ^ ~ t —^ ± nt- t —^-n% (3.5) 

Authors in Q and Q did not include the d k± terms in their calculations. We retain these terms as well and 
finally show that they contribute significantly to some of the form factors. All the individual distribution amplitudes 
are normalized to unity ie. J Q dx <f>\\ (x) = 1 etc. For the sake of completeness, we list the individual distribution 
amplitudes Q: 

4>\\(x) = 6.t(1 -x)[1 + 0.27(5[2a;- l] 2 - 1)] 
h<fi(x) = 3{2x - l) 2 + 0.3(2x - l) 2 [5(2a; - l) 2 - 3] + 0.21[3 - 30(2x - l) 2 + 35(2x - l) 4 ] 

dh^ s \x) 

h\l s \x) = — L — = Q(2 X - 1)[1 + 0.76(10x 2 - lOsc+l)] 

4> x (x) = 6x(l -x)[l + 0.3(5(2x - l) 2 - 1)] (3.6) 

q 

g { l\x) = -[1 + (2x - l) 2 ] + 0.24[3(2ir - l) 2 - 1] + 0.12[3 - 30(2x - l) 2 + 35(2x - l) 4 ] 
g { l\x) = Qx{\ - x)[l + 0.23(5[2x - l] 2 - 1)] 



„„, _ dg { l\x) 



dx 

Fourier transforming the wave functions with respect to the transverse momentum leads us to the wave functions 
expressed in the transverse separation 6-space. The variable b defines the transverse separation between the quark and 
the anti-quark inside the meson. In principle, the various DAs listed above can have different transverse momentum 
dependence. But for simplicity, we assume same dependence for all of them. This dependence is not known from 
first principles and therefore, it is generally assumed that the full wave function, including the transverse momentum 
dependence, is of the form 

*{x,k±) =4>(x)E(k 1 _) (3.7) 

such that apart from the normalization condition for 4> listed above, we have for the transverse part, J d 2 kj_E(fcj_) = 1. 
The functional form of is assumed to be a simple Gaussian distribution. In the &-space, we therefore have 

E P (6 P ) = exp - -e-M (3.8) 



For the heavy B-meson, following : 



and 0, 



we have the following momentum space projection operator 



*b(0 



75 



6 a 



where ^b(0 and ^b(0 are the B-meson wave functions corresponding to the DAs 4>b and (j)B defined as 



(3.9) 



(3.10) 



""2 ""2 

The term with coefficient A(Z) works out to be proportional to l + /tub and is therefore generally dropped. We 
retain this as well in the analysis and comment on its contribution later. Also, it has been argued that the \P 
term is sub-leading compared to the first term in the projection operator and is also dropped. This is one of the 
major issues of debate concerning the pQCD calculations. In the present study, we retain this sub-leading term also. 
Another reason of retaining this term is the confusion whether the pseudoscalar and axial matrix elements can both 
be described by the same DA in the heavy quark limit. It has been pointed out that the two DAs considered by the 
authors of \-\ (and of subsequent works based on their proposal), do not satisfy the equations of motion and further 
raise more doubts. In this study, we follow and use, as a model, the following form for the two DAs which satisfy 
the equations of motion 



2£ 2 m 2 



and 



m=\-— 



■ exp 



2 1 

exp 

7T LO B 



2u;% 



2lj% 



The quantity A(Z) introduced above is related to the DAs as follows 

A(Z)= / i+ d/(0 B -0+) 
Jo 

We next include the transverse momentum dependence for the B-meson through the following 

£b(&b) = exp 



uj%b 2 B 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



The question whether the above DAs are actually realistic and whether they correctly reproduce all the features is 
an intriguing one and has to be dealt separately. We therefore proceed with the assumption that they do or at least 
capture the important features. 



IV. SUDAKOV FORM FACTORS AND THRESHOLD RESUMMATION 



The resummation of large transverse separation regions leads to Sudakov form factors while the resummation over 
small fractional momenta leads to threshold resummation. 

As mentioned earlier, the Sudakov factor suppresses the double logarithms of the form In 2 Pb arising due to the 
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overlap of soft and collinear divergences. The parameter b regulates these kind of infra-red contributions. In terms 
of the variables 



q = In 



xQ 



V2A 



6= In 



the exponent of the Sudakov form factor, to NLO, reads 



A« q 



(q -b) + 



A (2) 



ln(2g) + 1 ln(26) + 1 



1-1 - 
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bK 



4/3? 



A« , / e 
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ln(2g) + 1 ln(26) + 1 



4/3? 

8/3? 

A( 2 )/3 2 
16/3f 



4(2)/32 ^_ t 

432|V [91n2(26)+61n(2 ' ) + 2] 



4^[ln 2 (2g)-ln 2 (26)] 



21n(2g) + 3 21n(2fe)+3 



A^(3 2 q - b 
16/34 S 2 



[2 ln(26) + 1] 



AWpj 
1728/9? 



181n 2 (2g) + 301n(2g) + 19 18 In 2 (26) + 30 ln(2b) + 19 

P 



where je is the Euler's constant and the various coefficients are given as follows: 



01 = 



33 - 2n/ 
12 



/3 2 = 



153- 19?i/ 
24 



(4.1) 



In [ % 



(4.2) 



(4.3) 



A^ 1 



,f2l 67 10 8_ , / 

A (2) = nt + -0i In 

9 3 27 ; 3 \ 2 



(4.4) 



(4.5) 



with rif being the number of flavours. The strong coupling constant to NLO is given by 

a,(fi) = 1 /3 2 lnln(/x 2 /A 2 ) 

it ftln^/A 2 ) /3? InV/A 2 ) 

The Sudakov factor therefore falls for large b regions and vanishes as b > 1/ A. Also, the wave functions (or equivalently 
DAs) defined above are valid only for scales 1/6. To have the correct results at an arbitrary /i a RG-evolution is 
required which generates an evolution factor in the exponential. We club this factor along with the Sudakov exponent 
and define the Sudakov as exp(-S) with 



S(x,b,Q,n) = 8(x,b,Q) + 8{l-x,b,Q)-^-hi: hll///A ' 



(4.6) 



01 ln(l/(6A)) 

In practical calculations, /i is identified with the factorization scale of the hard kernel. It is to be noted that for 
the B-meson, the Sudakov factor is considered only for the lighter quark. Therefore, the exponents of the individual 
Sudakov factors (including the evolution function) are given as follows: 



1 



S p = s(x, b p , m B ) + s(l - x, bp, m B ) - -5- In 



1 ln(l/(6 p A)) 



(4.7) 



10 

s -*^-4 ln w (48) 

The last major ingredient to be introduced is the threshold resummation. The double logarithms of the form In 2 x 
diverge at the endpoints and therefore they are also resummed. In order to achieve this goal, a jet function, St{x) is 
introduced which vanishes at the endpoints. It has been proposed 10] that for phcnomcnological studies, a simple 
parameterization can be used. The proposed parameterization is 

S t (x)= 2l+2 ^Ff / l + c) [x(l-x)r (4.9) 

with c ~ 0.3. We use this parameterization in our numerical study but agree that this may not be the accurate form. 
Finally, we remark that it is not very clear whether the power suppressed corrections arising due to working in 
different gauges for different quantities are going to play an important role [l^ . This is one feature that has to 
be checked but we leave it for a separate study. The DAs discussed in the previous section are assumed to be 
finally multiplied by the corresponding Sudakov factors, the transverse dependence carrying functions - the S's and 
threshold functions. 



V. EVALUATING B -> p FORM FACTORS 



The B — > p form factors as introduced above can also be written in a slightly different form. We make the following 
identifications. The coefficients of Pg + P p and Pg — P p define the form factors A 2 and A 3 , as discussed above. 
Instead, we define A2 and A2 as coefficients of Pb and P p such that A2 and A3 are simply sum and difference of the 
new form factors. The remaining two have the same definition as above. This has been done for convenience only. 
We evaluate the "parallel" and "transverse" contributions from both the diagrams separately. The individual ex- 
pressions read (again as mentioned in the last section, we assume the DAs to be multiplied by suitable Sudakov, 
threshold and momentum dependence factors and do not write them explicitly here in the expressions) 
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4\/2i77 2 a; 
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h\f\v)] 
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( 2iEf f 



"^V 2 t V^r/Vx = \/2 = 2r 2 x 

„2_^ *b(0 - mB(7? 2_ r 2) ^(0 ~ — *^(0 - — *b(0 



ra B {r\ 2 ~ r2 ) 



m B (r] 2 — 



*fl(0 + 



ni B (r] 2 ~ r2 ) 



*b(0 + 



4\/2i77 2 r 2 a; 
m B {rf — r 2 ) 



*b(0 



+ 2V2i/ ± / l ; s) (x) 



7= f 2 ^(0 - — *b(0 

^/2m B (i] -r 2 ) m B 



[x£rjm B - - fcx) 2 ]Nm| - fc 2 ] 
where the derivatives should be viewed as acting on the hard kernel and we have already set terms proportional to 
k± and l± appearing in the numerator to be zero. Similarly we have 

if B m 2 p m% 



M (2 l = (-4ttN c C f ) I liB ™P^ B J ( e * . n+ ) J dx #i ± d 2 k± as 
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*b(0 



1 77(77 + r 2 ) 
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+ 27/ ± /7; (s) (x) 
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These expressions are Fourier transformed to the 6-space. The wave functions/DAs are assumed to have been 
multiplied by the Sudakov, threshold resummation and momentum dependence factors. While Fourier transforming 
the expressions, we introduce the following hi functions: 
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and finally set r 2 = to retain terms up to twist-3. 

Before ending this section we would like to discuss the singular nature of the strong coupling constant as the scale 
A is approached. Close to this scale, it is easy to convince oneself that the perturbation theory should break down 
and therefore the calculation outlined above makes no sense as this scale is approached. For the reliability of any 
perturbative calculation, it has to be ensured that most of the contribution to the calculated quantity comes from 
the perturbative regime. It has been observed for the pion form factor that this is not really true and a large 
contribution does come from the non-perturbative region, something which makes the predictions highly unreliable. 
Near this scale, the coupling constant acquires a large value, thereby, rendering the breakdown of the RG improved 
perturbation theory. The infra-red behaviour of the strong coupling constant is one of the major challenges of modern 
day physics. However, it has been observed that there is a good phenomenological evidence to believe that close to 
(and below) the scale A, the coupling constant gets frozen to a value which is not a big number literally (see j^J] 
and references therein). The concept of a frozen coupling constant was proposed long ago by Cornwall |25|. Also, it 
has been proposed that the gluon propagator should be modified in order to circumvent the problem of singularity 
in the coupling constant and this has proved a successful phenomenological ansatz. If one adopts this point of view 
that as one approaches the scale A, the gluon propagator should be appropriately modified (using some regulator, 
say) and if one also assumes that the coupling constant gets frozen to some fixed value, then, one can try to carry 
out the analysis with the frozen value of the coupling constant for scales close to or below A. If the contributions 
are small, then this establishes the validity of the perturbative treatment. 



VI. BEHAVIOUR OF FORM FACTORS AND ROUGH NUMERICAL ESTIMATES 



The various form factors are read from the expressions obtained above. To numerically study the behaviour of the 
form factors, the choice of the scale /x has to be made. We take \i = Max{y/x^qmB 1 l/&s, 1/^p)- This scale must be 
greater than A in order to avoid the singular behaviour of the strong coupling constant. We would like to emphasize 
again that in this study we do not make any attempt of precisely calculating the form factors but focus on the 
behaviour of the same as functions of bs and b p and see whether the perturbative calculations make any sense. The 
choices for the other parameters are as follows: 

A = 0.25 CeF m B = 5.279 GeV m p — 0.77 GeV 



f B = 0.18 GeV f p = 0.198 GeV f± = 0.152 GeV 



lu b = 0.35 GeV uj p = 0.3 GeV 

The parameters uob and uo p should be 0(A) and therefore have been chosen as above. Figures 2 to 5 show the 
behaviour of various form factors for large recoil ie. q 2 — > or rj — > 1. From the figures it is quite evident that the 
Sudakov and threshold rcsummation factors are successful in regulating the endpoint behaviour of the form factors. 
Furthermore, it is not hard to convince oneself by looking at the figures that the perturbative region supplies the 
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dominant contribution. Also, it can be seen that for intermediate values of bs and 6 P , there is indeed an enhancement 
as can be expected and as discussed earlier. The Sudakov suppression is much weaker for the B-meson as is clear 
from the variation with bs ■ The reason for the appearance of a turning point while studying the variation with bs is 
the fact that the Sudakov suppression is weaker and there is a tendency for the strong coupling constant to blow up 
as 1/bs — > A. However, if one assumes the phenomenological ansatz of frozen coupling constant as A is approached, 
then explicit evaluations show that the contributions from the non-perturbative region are small. 
We summarize the numerical values for various form factors with these choices of parameters: 

V = 0.43 Ai = 0.39 



A 2 = 1.73 A 3 = 1.92 

Very evidently V and A\ are slightly overestimated while A 2 and A3 deviate significantly [2||. We would like to 
stress upon the fact that in this particular study we have not bothered to carry out a precise numerical analysis 
which would involve varying various other parameters like lob and uj p and then finally concluding which choices give 
the best results. It is straightforward to note that slightly higher values for both these parameters would have given 
more suppression. This has been the case with ^(J and [21! ] The best value which the authors in use for u>b 
was obtained by comparing the pQCD results to the pion form factor with the use of single wave function for the 
B-meson. However, when the other pieces are also included, this value is bound to change. 

In [2^, it has been remarked that the contribution from the A (I) term is roughly 20%. We agree with the general 
arguments and conclusions and reiterate that this term is crucial. Also, the ^b(0 term, though is sub- leading, but 
can still change the results when incorporated by at least a few percent. Furthermore, setting the two B-meson wave 
functions same is equivalent to setting the *&b(£) term to zero and also the A (7) term to zero. Therefore, the use 
of a single wave function for the B-meson does more than expected and therefore introduces more uncertainty than 
thought of. We attribute the reason of large deviations in A 2 and A3 to the presence of not really insignificant ^b(£,) 
and A(Z) terms. However, as mentioned earlier, the power suppressed terms that would arise due to different choices 
of gauges for defining different quantities have not been included and the relative sign of the same can be crucial in 
a complete and accurate numerical analysis. 



VII. CONCLUSIONS AND SUMMARY 



We have studied the B — > p form factors within the perturbative QCD approach including the twist-3 contributions 
to the DAs. For the B-meson, we employed a model for the wave functions that satisfies the relevant equations of 
motion and passes through this level of criticism. The use of single wave function has been criticised and we avoid this 
aspect by retaining the so called sub-leading pieces, which turn out to be not so sub-leading and insignificant. The 
appropriate Sudakov and threshold resummation functions have been included. With these ingredients the results 
in figures clearly are very encouraging and seem to suggest that the pQCD calculation of the various form factors 
is reliable and valid in the sense that dominant contribution clearly appears to stem from the perturbative regime. 
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However, since the aim was not a very accurate determination of the numerical values, we have overestimated, at 
least two of the form factors. But we believe that a careful choice of various parameters can lead to a consistent 
determination of the same. Also, there can be slight differences once a complete and precise jet function is employed 
in the calculations rather than the phenomenologically motivated function that has been used here. 
Encouraged by the results obtained, and mainly the fact that the endpoint behaviour is actually regulated and the 
contribution from the perturbative regime is the dominant one, we would like to propose the following as a working 
hypothesis for such calculations: 

* Use of two wave functions for the B-meson ie. retaining the generally thought sub-leading terms. 

* Including proper Sudakov and threshold resummation 

* To use, as a phenomcnological ansatz, the notion of a frozen coupling constant for scales very near or below 

The precise form of the B-meson wave functions is not known but we believe that the results will only improve with 
the advancement of our understanding of them. The results of this study, particularly the behaviour of the form 
factors with bs and b p are very encouraging and require a detailed and careful numerical scrutiny. 
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FIG. 4: Variation of A 2 with b B (left) and b p (right) 




FIG. 5: Variation of A3 with bs (left) and b p (right) 



